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Abstract 

The Landau-Yang theorem, forbidding transition amplitudes be¬ 
tween a massive spin-1 particle and two photons, is widely assumed 
to apply to other massless spin-1 hnal state particles as well. We 
show that this is not true in Standard Model QCD, so that for in¬ 
stance antisymmetric colour-octet spin-1 quarkonia can be formed 
by two on-shell gluons at O^a"^) in perturbative QCD. 
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1 Introduction 


The Landau-Yang theorem mm is an old and well-established result in QED, 
forbidding the decay of a massive spin-1 particle into two photons. A fine exam¬ 
ple is provided by orthopositronium: since it cannot decay into two photons but 
only into at least three, it has a much longer lifetime than spin-0 parapositro- 
nium. Over time the Landau-Yang theorem has come to be interpreted more 
broadly to imply that the decay of a spin-1 particle into any two massless spin-1 
particles is also forbidden, and this is used in studies of quarkonium production, 
such as for instance in [3] under the name of the generalized Landau-Yang selec¬ 
tion rule. In this article we present a very simple derivation of the Landau-Yang 
theorem based on Bose statistics, point out a loophole presented to non-Abelian 
gauge theories and prove that the generalized Landau-Yang theorem does not 
hold for at least some QCD processes. 

We give explicit results for next-to-leading order (NLO) processes and explain 
the reason behind the fact that a generalized Landau-Yang theorem seems to 
exist at leading order (LO), but is violated at NLO, as hrst noted in |3]. Our 
results are general and we show that they reduce to the approximation given in 
the aforementioned paper when the appropriate limits are taken. 


2 The Landau-Yang theorem in QED 

Before turning our attention to the validity of an extension of the Landau-Yang 
theorem to the case of QCD, we first revisit the original theorem due to Landau 
and Yang within the framework of QED. 

The method Yang used to obtain his result is based on the transformation prop¬ 
erties of the annihilation and creation operators of the electromagnetic field un¬ 
der rotations and parity transformations. These lead to a set of selection rules 
which rule out the spin-1 decay into two photons. We present here a different 
method using the language of transition amplitudes. 

Consider a spin-1 particle of mass M, momentum P and polarization eg de¬ 
caying into two photons with momenta qi and q 2 and polarizations ei and £2 
respectively. This means that we have P • eg = gi ■ Ci = (?2 • £2 = 0. We also 
define 



92), 


( 1 ) 


so that P ■ r = 0. Hence P and r are independent, orthogonal momenta. 

Owing to the Ward identity we are allowed to add to any polarization vector a 
term proportional to its corresponding momentum. To be specihc, we re-gauge 
£ 1^2 as follows: 


11 _ n 

m = ef 


(e» ■ qk) 

{qt ■ Qk) 


9 * 


for (j, fc) = (1,2) or (2,1). 


( 2 ) 
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These have the additional property P ■ 771^2 = f ■ Vi,2 = 0 . By simple exhaustion 
it is easy to see that the decay amplitude ^ must necessarily be of the form 

^ = ai{r ■ eo)ir]i ■ 772) + 02 e(f’, eo, 771,772) + asir ■ eo)e(P, r, 771,772), ( 3 ) 

where e{p,q,k,l) represents the contraction of its arguments with the Levi- 
Civita symbol: e{p,q,k,l) = e^'^^’^p^qi^kpla- The coefficients ai can apart from 
constants only depend on M. 

The Landau-Yang result is now easily obtained by interchanging the two outgo¬ 
ing photons (771 -O' 772,r —>■ —r). Under this transformation each term changes 
sign, while their coefficients do not. To obey Bose symmetry, the coefhcients 
Qi all have to be equal to zero, which is the Landau-Yang result. 

The loophole in this argument for the case of QCD is of course that the coef¬ 
ficients Gi contain a colour structure as well. The colour-antisymmetric part of 
these coefficients therefore leads to an amplitude that is Bose-symmetric and 
hence not forbidden. In many cases, such as the decay Z ^ gg involving quark 
triangle diagrams, the amplitude is symmetric in the gluon colours and the gen¬ 
eralized Landau-Yang theorem holds. This has been known for a considerable 
time, see for instance [S] where theories beyond the Standard Model are con¬ 
sidered in which colour-antisymmetric effective interactions can be constructed. 
The aim of this article however is to show that also in the minimal Standard 
Model the generalized Landau-Yang theorem can be evaded in QCD and to give 
an explanation of how this works (both at LO and at NLO). 


3 The Landau-Yang theorem in QCD at LO 

We now turn to the case of QCD by replacing the photons in the previous section 
by gluons. As stated before, the coefficients now contain colour structures 
which might provide a loophole to evade the Landau-Yang theorem in QCD. 
We consider the process 

5(92,772, fc) qipi,a) q{p2,b), ( 4 ) 


where two incoming gluons with momenta 771^2, polarization vectors 771^2 and 
colours j, k produce a quark-antiquark pair with momenta pi^2 and colours a, b. 
The outgoing quarks have mass m. At the kinematic threshold this process is 
equivalent to the decay of quarkonium into two gluons considered in the rest- 
frame. 

At LO in the strong coupling gs, this involves three diagrams as depicted in 


Figure 1 In order to discuss the amplitude, we introduce 


P= ^(91 +92) = ^{Pi +P2), 


9 = ^{Pi-P2), 


r = ^{qi - 92)- 


By also defining the slightly altered Mandelstam variables 

s = 2(gi • ( 72 ) = V, ^ =-2((7i-P i), u= -2{qi-p2), 


( 5 ) 


( 6 ) 
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Figure 1 : The diagrams contributing to gg —)■ qq at LO. 


we can express the three diagrams referred to as the s-, t- and u-channel graphs, 
as 


^o,s = - ^u{pi){r]i ■ r] 2 )fv{p 2 ) X C, ( 7 ) 

. 2 

^o,t = - ■ m)% + ( 9 -^2)1^1 -^if%)viP 2 ) X {A + C), (8) 

. 2 

•^o,u =-^-^u{pi){{q-p 2 )il^ + {q-gi)il^ + TllJil^v{p 2 ) X {A-C), ( 9 ) 

where C = r^]“^ and ^ = {T^, denote the commutator and anti¬ 

commutator of the colour generators and in the fundamental representa¬ 
tion. Adding them yields the total amplitude at lowest order 

-^0 = ^r—.u{Pi) r v{p2) y. [a + ^— -C), 

2 ut \ s y 

r = (g • ?7 i)i^ 2 + - 4fy£(p,?7i,r,772)7® -b ■ V2)f- (lO) 

Note that the part proportional to the commutator of the colour matrices van¬ 
ishes at the kinematic threshold (where t —>■ u and g —> 0) because of its prefactor 
(t — u) / s. This means that the amplitude becomes colour symmetric and there¬ 
fore completely analogous to the case of QED. 

The only surviving term at the kinematic threshold is the one involving 7^. This 
structure, however, vanishes when we consider the outgoing quarks to be in the 
spin -1 ortho-quarkonium state, which leads to an extension of the Landau-Yang 
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theorem to QCD at lowest order. 


The specific factorization structure exhibited by Eq.(lO), where the symmetric 
and antisymmetric colour part differ by a factor (t — u)/s, arises due to a cancel¬ 
lation between the s-channel and the t- and u-channels. It is this factorization 
structure that is responsible for the vanishing of the colour antisymmetric part 
of the amplitude. 


Before performing a NLO calculation for the process in which a quark-antiquark 
pair is produced by two gluons, we list all possible spinor structures and cate¬ 
gorize them according to their parity, their survival in case the two quarks form 
a spin-1 ortho-quarkonium state and whether they survive at the kinematic 
threshold or not. 


Structure 

Parity 

Threshold 

Spin -1 

u{pi)v{p 2 ) ivi ■ m) 

-b 

- 


u{pi)v{p2) ivi ■ q){m ■ q) 

-b 

- 


u{pi)j^v{p 2 ) im ■ m) 

- 

-b 


u{pih^v{p2) im ■ q){f]2 ■ q) 

- 

- 


u{pi)il)^v(jp2) im ■ q) 

-b 

- 


uiPi)fl^v{p2) im ■ q) 

-b 

- 


u{pi)fv{p 2 ) im ■ m) 

-b 

-b 

-b 

u{pi)fv{p2) im ■ q){v2 ■ q) 

-b 

- 


u{pi)i^'i/i^v{p2) im ■ q) 

- 

- 


u{pih^fi2'>j{p2) im ■ q) 

- 

- 


u{pi)i^fv{p 2 ) ivi ■ m) 

- 

-b 


u{pi)'y^fv{p2) im ■ q}(g2 ■ q) 

- 

- 


uipi)$Jvip2) im ■ q) 

-b 

- 


uipi)^^fvip2) im ■ q) 

-b 

- 


uiPlh^'^^fvip2) ip2 ■ q) 

- 

- 


uipi)-i^'rii^fvip2) im ■ q) 

- 

- 


uiPliiTfl^T/l^ - Tfl 2 lfi^)'f^iP 2 ) 

-b 

- 


uiPlh^i-^J^ - 'j 2 ^Jvip 2 ) 

- 

-b 



-b 

-b 

- 


Table 1 : Categorization of possible spinor structures in gg —> qq 


Every spinor structure can be reduced to the ones found in [Table 1 | by anti¬ 
commuting the various vectors or using identities like 


■ ^2)/ - (r ■ C m,‘i 72 ) 7 d, (11) 
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where E stands for the center-of-mass energy of any one of the particles. We 
have used the Pauli-identity and the fact that p is the only parameter with a 
time-component and no space-components. 

A plus sign in [Table l] means that the spinor structure exhibits the desired 
behaviour. For parity, a plus sign means that the structure is even under a 
parity transformation. For the threshold entry a plus sign indicates a structure 
that does not vanish at the kinematic threshold. A plus sign in the final entry 
indicates that the specific structure does not vanish when the two quarks are 
considered to form a spin -1 ortho-quarkonium state. 

From [Table l| it appears that only one spinor structure is possible for our pro¬ 
cess of interest if the quarks are to form an ortho-quarkonium state at the 
kinematic threshold. This structure, u{pi)fv{p2){rii • 772), is also present in the 
LO amplitude but vanishes due to a prefactor {t — u)/s. 

4 Gluon-initiated quark-antiquark pair produc¬ 
tion amplitude at NLO 

The last term in the Born amplitude (Eq. proportional to {rji ■ 772)/, only 
vanishes at the kinematic threshold due to its prefactors, while the other terms 
vanish by themselves. This means that there might be a non-vanishing transi¬ 
tion amplitude involving this term at higher orders. 

To monitor divergences at the kinematic threshold, we introduce the threshold 
parameter 


I 4 777^ 

which approaches zero at the kinematic threshold. All spinor structures are ex¬ 
panded in terms of this parameter, as are the Mandelstam variables and other 
prefactors. 


To obtain the NLO amplitude for quark-antiquark pair production by two glu¬ 
ons, all graphs are considered at the kinematic threshold, except for the Coulomb 
graphs given in [Figure 2 Diagrams containing infrared divergences, in which a 
soft gluon is radiated by one of the external particles, do not need to be taken 
into account because they can be shown to be proportional to the (vanishing) 
LO amplitude. The expression for each graph is then reduced to scalar integrals 
by the method of Passarino-Veltman reduction mzii. The Coulomb graphs 
give rise to scalar Coulomb integrals which are proportional to 1//3 and need to 
be expanded in terms of the threshold parameter. After all 1//3 divergences are 
cancelled by terms arising from the expansions of the prefactors (containing q, 
s, t and u) of the Coulomb integrals, the kinematic threshold limit can be taken 
for the Coulomb integrals as well. When all contributions are added, the fact 
that we are considering the quark-antiquark pair to be in the ortho-quarkonium 
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Figure 2: The Coulomb diagrams for gg —>■ qq. 


state is taken into account and we obtain the final resullQ 


9tF 

192TO7r 


— 2 ln(2^) H- \- in 




X (^\g2\cos{(j)2)u{pi)'ilj^v{p2) + \-ni\cos{(t>i)u{pi)Tli^v{p2) 

gtF 


2 ln(2^)-2 — in 


192m^7r 


cos(6») X u{pi)fv{p 2 ){gi ■ 92 ) 


^9s 


\92w?n'^ 

i 


^ Quarks 

+ F( 14-71^-28In 


E . / . ,110^ 111'^ »1 o ,110^ 

F' 6/3' ln(a:'J^ + 12^ + 2 + - 

X u(pi)/7;(p2)(»7i ■ g2 )[ T \ T ^ Yf ^. (13) 


Here we have defined e by d = 4 — 2e, where d is the number of space-time 
dimensions. Furthermore we have introduced 


/3' = 1 

/ 4m2 

V s + zd’ 

= 

F= I 

f TO^e'’'® ^ 

and 

dTT/r^ y 


/3' + l’ 




47r/i^ 


( 14 ) 


“^The complete calculation can be found in [^. 
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where rUg denotes the mass of a quark occurring in a quark loop and /i indicates 
the ‘t Hooft scale. The angles 4'i,2 refer to the angles between pi and 771^2 and 
the angle 6 denotes the angle between the incoming and the outgoing particles. 
Strictly speaking these are not defined at the threshold; this is a reflection of 
the fact that a realistic treatment would involve folding the amplitude with a 
quarkonium wavefunction. Secondly the colour symmetric part of the amplitude 
does not vanish, which might seem to contradict the Landau-Yang theorem for 
QED. However, from the very fact that the angles 4>i^2 and 9 occur, it can 
be seen that these contributions come from partial waves beyond the S-wave. 
Therefore the Landau-Yang theorem should not be expected to hold for these 
terms. In fact, the colour symmetric structures in Eq.(13l arise as a result 


of Coulomb corrections proportional to 1//3 that pick out /3-suppressed terms. 
These Coulomb corrections are formally part of the qq bound state, because 
they appear in the limit where the slow moving, outgoing quarks exchange a 
virtual gluon. Because of this, only the last part of the NLO result describes 
the amplitude of creating a fundamental spin-1 particle, while the other terms 
belong to the description of the formation of a bound state from two individual 
quarks. As the Landau-Yang theorem explicitly deals with fundamental spin-1 
particles, we focus only on the colour antisymmetric part: 






Quarks 


192m27r2 

+ -28\n{2)j 


In o ' 


X u{pi)fv{p 2 )i'ni ■ m)[T\T'^] J,, (15) 


in the following discussion. The F and F' in this expression can in principle be 
set to 1, since there are no 1/e poles present in 


5 On the fate of the Landau-Yang theorem in 
QCD 

The term in the Born amplitude that might have provided a loophole for the 
Landau-Yang theorem in QCD is given by 

*^0 = —^m(pi)/w(p2)(??i • 97 ) X ^— - Co - (16) 

2 ut s s 

As {t — u)/s ^ /3, this term is doubly suppressed at the kinematic threshold. 
The reason for this suppression, as stated before, is the factorization due to a 
cancellation between the s-channel and the t- and u-channel graphs. 

We now consider the total amplitude squared up to order 

F F F F F .^ 2^0 . ( 17 ) 


where the subscript denotes the order. 

The NLO amplitude has no means of cancelling the suppressing terms of the 











Figure 3: The NLO graphs involving quark loops. 


Born amplitude, which means that the first three terms in Eq.(17) drop out 


at the kinematic threshold. Concerning the fourth term, we will only focus on 
the colour antisymmetric part -# 5 ^ as discussed in the previous section. The 
part that involves a sum over quarks will lead to a double sum over quarks in 
In order for +^ 2 *^ 0 *) cancel these terms, next-to-next- 

to-leading order (NNLO) diagrams involving two quark loops are required. Such 
NNLO diagrams do not give rise to integrals containing factors of 1//3 which 
could cancel the suppressing factors of This means that the terms propor¬ 
tional to F' in the NLO amplitude will be present in and the Landau-Yang 
theorem does not hold in QCD! 


The terms proportional to F' originate from the graphs involving quark loops, 
given in Figure 3 These three sets of graphs are all s-channel diagrams, which 


means that the cancellation between the s-channel and the t- and u-channel 
that gave rise to the specific factorization structure at LO does not occur for 
the graphs given in |Figure~3l Thus the loophole for the Landau-Yang theorem 
in QCD, which was closed at LO by threshold suppression, opens up at NLO. 


To summarize, the calculation of the 1-loop amplitude for quark-antiquark pair 
production by two gluons involves, among others, terms proportional to the sum 
over virtual quarks. We argue that, at the cross-section level, these terms can 
not be cancelled by terms resulting from the interference between the Born am¬ 
plitude and the two-loop amplitude. As all lower order contributions vanish, we 
conclude that for this process, the generalized Landau-Yang theorem is violated 
at two-loop level in the cross-section. 

We thus have here a mechanism to form vector quarkonia in an antisymmetric 
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colour-octet state using on-shell gluons, which can be described by hadronic 
PDF’s, without resorting to input from outside the Standard Model. 

Once a particular flavour for the outgoing quark-antiquark pair is chosen, quite 
often an effective-field-theory (EFT) motivated approximation is employed. Heavy 
quarks, i.e. quarks with a mass larger than the mass of the produced quarks, are 
treated as decoupling. Therefore, the sum in only involves quark masses 

rriq up to the mass of the produced quarks in that approximation. 

The light quarks with a mass smaller than the mass of the produced quarks 
are treated as massless, causing all terms inside the sum to vanish except for 
the constant ‘2’. For the quark flavour equal to the flavour of the produced 
quarks, we have = to, which leads to a factor of 14 — 37r^/2. In this EFT 
approximation we therefore have 

+ 2-711 + (lA - TT^ - 281n(2)) 

X u{pi)fv{p2){vi ■ m)[T^ (18) 

where ni denotes the number of light quarks. 

The non-vanishing NLO amplitude for the process of charm-anticharm produc¬ 
tion by two gluons and the fact that the Landau-Yang theorem does not hold in 
this case was first noted in |3]. Our result (Eq.([I^) represents a more general 
process in which any quark-antiquark pair can be produced by two gluons and 
where the contributions from individual quark loops are not subject to approx¬ 
imations. 

When we consider the outgoing quark-antiquark pair to be a charm-anticharm 
pair (to = TOc), our expression for the NLO amplitude indeed reduces to the 
result obtained in [1] if we treat the lighter quarks as massless and the heavier 
quarks as decoupling. 


= _ iRs _ 

NLO.approx 192^2,^2 
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